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Abstract 

We show that a surface term should be added to the Einstein-Hilbert action in 
order to properly describe quantum transitions occurring around a black hole. The 
introduction of this boundary term has been advocated by Teitelboim and collaborators 
and it has been used in the computation of the black hole entropy. Here, we use it to 
compute the gravitational corrections to the transition amplitudes giving rise to Hawking 
radiation. This surface term implies that the probability to emit a particle is given by 
e -AA/4 w h ere /\^4 j s change in the area of the black hole horizon induced by the 
emission. Its inclusion at the level of the amplitudes therefore relates quantum black 
hole radiation to the first law of black hole dynamics. In both cases indeed, the term 
expressing the change in area directly results from the same boundary term introduced 
for the same reason: to obtain a well defined action principle. 
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1 Introduction 



There are two possible approaches to the gravitational back reaction to Hawking radi- 
ation. The first is to develop a microscopic theory of quantum gravity and to use it 
to calculate the properties of black holes. This program has been partially realized in 
the context of super-string theory [0. The second is to use Hawking's calculation as a 
starting point to compute gravitational corrections to black hole evaporation. Hopefully 
these two approaches should meet in some middle ground. 

In the more conservative second approach, the back reaction has been addressed along 
two complementary lines. The first is the "semi-classical" theory wherein the metric 
remains classical, but one takes as source in Einstein's equations the mean (average) 
energy-momentum of the second quantized fields propagating on this self consistent 
background 0]. This describes the mean back reaction to Hawking radiation and 
confirms that black holes evaporate adiabatically in a time 0(M 3 ) as first predicted by 
Hawking 0. 

The second line of attack is to take into account the dynamics of gravity at the level 
of transition amplitudes before performing any quantum averaging. In the approach 
initiated by Keski-Vakkuri, Kraus and Wilczek (KKW) ||, one replaces the matter 
action S mat ter in a given geometry by the self consistent action of matter plus gravity 
Smatter+gravity One then postulates that the wave functions governing transition am- 
plitudes are WKB, i.e. are of the form e lSm+g . This approach closely corresponds to 
the computation of the transition amplitudes in quantum cosmology performed in ||. 
Indeed, the gravitational waves appearing in matrix elements are WKB solutions of the 
Wheeler-De Witt equation and thus have also their phases determined by the matter 
energy of the quanta involved in the transition. Another useful analogy to this approach 
is provided by the quantum description of the trajectory of an accelerated detector ex- 
periencing the Unruh effect 0. In that case, the given trajectory (which plays the role 
of the classical geometry in black hole physics) has been replaced by WKB waves so as 
to take into account recoil effects according to Feynman rules@||. 

The aim of the present work is to clarify the mechanisms at work in the approach of 
KKW. In particular, we determine what is the correct gravitational action that must be 
used upon describing black hole radiance. There is indeed some ambiguity in the choice 
of the boundary terms at the horizon. On physical grounds, we select the appropriate 
one. This amounts to work with the boost parameter defined at the horizon as the time 
parameter instead of the time defined at infinity. In this, our treatment of the boundary 
terms closely follows the work of Banados, Teitelboim and Zanelli|TO|. 



In a former paper |IT|, we already used their improved action in a Euclidean context. 
Our aim was to compute self consistently the action of instantons governing probability 
transitions in the presence of horizons. The main advantage of the Euclidean approach 
lies in its simplicity. In that case indeed, the choice of the correct action is particularly 
transparent. However this approach has an important weakness: one must postulate, 
as in all works based on the Euclidean path integral see e.g. [O, that no conical 
singularity be present in the Euclidean geometry. This condition should be proven 
from first (quantum) principles and not assumed. Remember that in the absence of 
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gravitational back-reaction it is the fact that the quantum state is vacuum that implies 
Green functions possessing well defined Euclidean properties. 

The present work, that of KKW, and that of JTT| all reach the same conclusion. By 
using the Einstein-Hilbert action rather than only the matter action as the phase of 
WKB waves, one finds that the rate of emission of particles of energy A from a black 
hole of mass M is 

R M ^ M _ x = N(\,M)e- AA ^ M ^ (1) 

where AA(X, M) = A(M) — A(M — A) is the change of area of the black hole due to 
the emission of a particle of energy A, and N is a phase space (also called grey body) 
factor which cannot be calculated in this approximation scheme. To first order in A, 
AA(\, M)/4 = X/8nM and one recovers Hawking's canonical result characterized by a 
temperature. 

One of the main interests of eq. ([T]) is to suggest that the thermal properties of 
Hawking radiation and the first law of black hole mechanics both stem from the same 
fundamental principle valid at the quantum level, i.e. before averaging. This would 
explain why they give a consistent description of the thermal properties of black holes. 
We recall that the first law relates neighboring black hole classical configurations and 
as such seems to be completely independent of quantum black hole radiance. A first 
direct indication of the relationship between the first law and black hole evaporation is 
provided by the instantonic approach ||1 1||. Indeed in this work, AA(X, M)/4 has exactly 
the same origin as the A/ 4 term in the partition function of a black hole calculated using 
Euclidean path integral techniques ||13|| |TT| . In both these cases, this term arises from a 
boundary term which is introduced at the horizon in order to select the Hamiltonian 
variational principle appropriate for the physical problem under investigation. And in 
the first law as formulated by Wald |TJ| , dA /4 is also equal to the variation of the same 
boundary term at the horizon. At this point it should be noticed that this boundary term 
does not appear in the calculation of KKW. Our aim is therefore to carefully reformulate 
the calculation of transition amplitudes to establish that the boundary term must be 
present right from the beginning. 

We finally note that the present procedure applies to all processes occurring in the 
presence of horizons and leads to transition rates determined by the corresponding 
changes in area. In particular it also applies to the transitions of accelerated systems. 
This suggests that there might exist a quantum statistical extension of the thermodynam- 
ical approach to general relativity [|1J which is valid everywhere and not only applicable 
to black hole horizons. In this case one would obtain a kind of universal "holographic 



principle" []T6| valid for all horizons and expressing the changes of the effective degrees 
of freedom to which systems living on one side are coupled. 

The plan of the paper is as follows: we begin by reviewing the standard derivations 
of Hawking radiation without back reaction, following || and 0. We then review the 
role of boundary terms in the Einstein-Hilbert action. Finally we rederive eq. ([!]) in two 
complementary ways. The first emphasizes the role of detectors which measure the pres- 
ence of particles and which deform differently the geometry according to their quantum 
state. The second way only makes use of the regular structure of the wave functions 
near the horizon when expressed in advanced Eddington-Finkelstein coordinates. 
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2 Hawking Radiation 



In this section we fix the notations and review two standard derivations of black hole 
radiation. In the first, following Unruh's workf^], one introduces a particle detector 
whose position is fixed and one determines the populations of quanta from its transition 
rates. In this way, one only uses basic quantum mechanical rules. The second approach 
is more intrinsic and closer to the original Hawking's derivation @. Black hole radiation 
is established through the Bogoljubov transformation relating in-modes which determine 
the Heisenberg state and out-modes which define the on-shell particles found at infinity. 
The specification of the in-modes will be obtained by imposing well defined analytical 
properties on the horizon[]T|][|T7[ . 



The metric of a Schwarzschild black hole is powerful 

ds 2 = -(1 - ^)dt 2 + (1 - —^dr 2 + r 2 dQ 2 (2) 

We introduce the light like coordinates u and v: 

v,u = t±r* , r* = r + 2M\\n(r - 2M)\ (3) 
and the Kruskal coordinates U and V: 

U = --e~ KU , V = -e KV (4) 

K K 

where k = 1/4M is the surface gravity. 

For a black hole formed by the collapse of a star, the outgoing modes solution of the 
Klein Gordon equation have the following form near the horizon: 

e~^ u Y lm (Sl) 

0w,l,m = r . (5) 

V47rtJ r 

Further from the horizon this is not an exact solution because of the potential barrier 
which surrounds the black hole. For simplicity in this article we shall neglect the poten- 
tial barrier. It would only encumber the expressions, whereas all the physics lies near 
the horizon where it plays no role. 

A way to understand the structure of the modes eq. (||) is that they correspond 
to the absence of particles as seen by an infalling observer. Indeed the proper time 
lapses of an infalling observer are proportional to AU. This is directly apparent by re- 
expressing the Schwarzschild metric in U, V coordinates. For r — 2M 2M, one finds 
ds 2 ~ —dU dV + r 2 d 2 Q which explicitizes the initial character of U near r = 2M. 

The field operator can be decomposed as: 

$ = a u,i,m4>u,i,m + h.c. + ingoing modes (6) 

ui,l,m 

By definition, the Unruh vacuum state \0u > is annihilated by the a w ; jm operators. 
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Consider now a particle detector at fixed radius r and angular position f2 which has 
two levels |e > and \g > of energy E e and E g with AE = E e — E g > 0. Its coupling to 
the field $ is given in interaction representation by 

H mt = j${t,r,n) (e +tAEt \e >< g\+h.c.) (7) 

where 7 is a coupling constant. Initially the detector is in its ground state, hence the 
combined state of the detector and field is \0jj > \g >. In the interacting picture, at late 
times and to first order in 7, the state of the detector plus field is: 

\i/j g > = \Q V > \ g > _i 7 f dte iAEt 3>(r,n,t)\0u > |e > 

= 10* > \g > -vy E / dte iAEt e-^- Kt ^&at lm \0u > \e > (8) 



ojlm 



where C = (l//«)exp(Kr + |ln(r — 2M)). Similarly if the detector is initially in its 
excited state, the amplitude of the detector plus field at late times is 

>= |CV > |e > -i 7 E / cite-^e-^-^^Sa^lOc; > |p > (9) 



ujlm 



The ratio of the amplitude of getting excited to the amplitude of getting deexcited is 

A , f rlfp +iAEt P -iuCe~ Kt 

~A ] _ f rlfp~i/\Et p ~iujCe- Kt \ ' 

•ST-e—>g,u),l,m J o 

By replacing t by t+in/n in the upper integral, one obtains that A g ^ e = Al^ g e~ 7TAE ^ K , 



sec 



|I8~| . Thus the ratio of the probabilities of transition is: 



\<A-e— >g,uj,l,m\ 



corresponding to the rates in a thermal bath at temperature k/2tt = l/8nM. In equi- 
librium the ratio of rates is equal to the ratio of probabilities to be in the excited or 
ground state. This derivation is equivalent to calculating the Bogoljubov transforma- 
tion between Unruh modes and Schwarzschild modes. Moreover, it provides a physical 
interpretation to individual Bogoljubov coefficients as transition amplitudes ]7| 

In the second way, the spectrum of emitted particles is determined in terms of modes 
which are eigenvectors of id% = A. There are two modes for each value of A which are 
non vanishing either outside the horizon or inside the horizon: 

¥W = e ^ m( "W -2M)) (12) 

V 47TA T 

Where we have once more neglected the potential barrier outside the black hole. The 
Schwarzschild modes tp + are those which would naturally be used by a static observer 
around the black hole to describe the presence or absence of particles. On the other 
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hand the modes 0A,z,m associated to the Unruh vacuum state are linear combinations of 
ip + and <£j_. To determine this linear combination one must impose that the 0A,z,m have 
only positive frequency for an infalling observer. One way to construct the <p\ modes 
is to express the Schwarzschild modes in Kruskal coordinates: <^a± ~ (^fU) lX ^ K . Since 
AU is proportional to the proper time of an infalling observer, <fi\ must be the linear 
combination of tp\± which is analytic and bounded in the lower half of the complex U 
plane: 

where the normalization factor ensures that the </> modes have unit norm. The weights 
in eq. ( P^) define the Bogoljubov coefficients ot\ and (3\. Their ratio satisfies 

This implies that the Unruh vacuum is a thermal distribution of Schwarzschild particles 



at temperature k/2ttM (compare with eq. fllTD ) 



We will find it convenient below to use the same argument, but rephrased in Eddington- 
Finkelstein coordinates v, r, Q in which the metric has the form: 

ds 2 = -(1 - 2M/r)dv 2 + Idvdr + rW (15) 

Near the horizon the metric takes the simple form ds 2 ~ Idvdr + r 2 dQ 2 which shows 
that v, r are inertial coordinates near the horizon. An infalling observer follows, near the 
horizon, approximately the trajectory v = const, so that — r is approximately a proper 
time parameter along his world line. <f)\ t i tm must have only positive frequency in — r at 
fixed v . Near the horizon, in Eddington-Finkelstein coordinates, the (p modes have the 
form (fx,i, m ,± — e~ l ~( KV ~ ln ( r ~ 2M ^ 9 (±(r — 2M)). The regularity condition imposes that 
to obtain a mode, one analytically continues tp + in the upper half complex r plane to 
obtain the relative weight of the ip + and </?_ modes, as shown in [17] . In this way one 
recovers eq. (O). 



3 Boundary terms in the Einstein-Hilbert Action 

In this section we examine the boundary terms in the Einstein-Hilbert action. We start 
our analysis with the action expressed in Hamiltonian form: 

S = Jdt {ir ab g ab +P q-NH- N l H t } (16) 

where g a b is the spatial metric, n ab its conjugate momentum, q and p the coordinates 
and momentum of matter, iV and N 1 the lapse and shift, and H and Hi the energy and 
momentum constraints. 

We shall consider this action defined in the right quadrant of the full Kruskal mani- 
fold, and take the time foliation to resemble Schwarzschild time. The equal time slices 
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then have two boundaries: one at spatial infinity and one at the horizon (more precisely 
at the intersection of the past and future horizon). Variations of the action eq. ([T6| ) gives 
rise to boundary terms because of the presence of spatial derivatives of the metric in H 
and Hi. In order that the action be extremal on the solutions of the equations of motion, 
one must add surface terms to S that cancel the boundary terms. We now review these 
boundary terms, starting with the boundary at infinity. Details of the calculations will 
not be presented. They can be found in in many papers, see for instance [|TI| [^DJ []ID| 1 



Varying eq. ( |lq) while imposing that the space time is asymptotically flat yields 

SS = J dt terms giving eqt. of motion + t^SMADM (17) 

where too = / dtN is the proper time at infinity and Madm is the mass at infinity 
(defined in terms of the behavior of g rr for large r). In this form the action is extremal 
on the equations of motion if one varies among the class of metrics for which Madm is 
kept fixed. 

On the contrary, if one varies among the class of metrics for which too is kept fixed 
but Madm is arbitrary one must subtract the surface term Madm^oo to the action S so 
that the variation of the new action 5" yields 

SS' = J dt terms giving eqt. of motion — Madm^oo (18) 

This second form of the action is generally more convenient because one can calculate S' 
as a function of too, and then use this expression to calculate the time evolution of the 
matter and metric. Indeed since Madm is a conserved quantity, dS' / 8Madm = const 
gives the time evolution. A second more physical reason is that one may want to compare 
the evolution of systems with different ADM masses, and hence one does not want to fix 
it to start with. However we must emphasize that there is nothing fundamentally wrong 
with the variational principle based on S. This will be important below. 

We now turn to the boundary term at the horizon. The analysis proceeds exactly 
in parallel with the proceeding one. The form of the boundary term is dictated by the 
fact one requires that near the horizon the lapse and shift vanish (or equivalently that 
at p = 0, the momentum ir? and the derivative of the area of the surfaces of constant p, 
d p A, vanish, see pl| ) and hence the metric can be put in the form 

ds 2 = -p 2 n 2 dt 2 + dp 2 + r 2 (p)dn 2 (19) 



For a Schwarzschild black hole p = ^8M(r - 2M), k = 1 /AM and t = t c 



Upon varying eq. (16) one then finds a boundary term at the horizon 



SS = J dt terms giving eqt. of motion — QSA/8n (20) 

where = J ndt is the hyperbolic angle and A is the area of the horizon. The action in 
this form is therefore extremal on the equations of motion for the class of metrics which 
have fixed horizon area. 
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One can add to this action a surface term at the horizon and define S" = S -\- QA/8tt 
so that its variation takes the form: 

5S" = J dt terms giving eqt. of motion + A5Q/8n (21) 

This action is extremal on the equations of motion provided one makes variations in the 
class of metrics which have fixed at the horizon. 

This is the action that must be used when calculating the action of the Euclidean 



continuation of the black hole|T0[- Indeed regularity of the Euclidean manifold at the 
horizon imposes that the Euclidean angle = 2tt. The surface term in the action then 
contributes a term A/4 to the partition function which is interpreted as the entropy of 



the black hole. This is also the action that was used in |TT|| to calculate the self-consistent 
actions of instantons. In the present paper, we shall use it in the Lorentzian sector. The 
new aspect brought in by its use is that the time evolution will be given in terms of the 
boost parameter B rather than the time at infinity. 

Thus upon considering processes occurring around a black hole, there are a priori 4 
actions that can be considered according to choice of the surface terms. Two however 
are unphysical. Indeed fixing both the ADM mass and the horizon area is inconsistent. 
For instance for the vacuum spherically symmetric solutions, fixing the ADM mass de- 
termines the horizon area. Similarly fixing both the time at infinity and the hyperbolic 
angle O is inconsistent. Thus one is left with two possibilities: fixing the ADM mass 
and G or fixing and A. A more mathematical reason why these are the only two 
possibilities is that the constraints H = and Hi = viewed as differential equations 
need boundary conditions in order to yield a unique solution and fixing the ADM mass 



or the horizon area but not both provides the required boundary data|^T • 

The choice among these two possibilities is dictated by physical considerations. If the 
ADM mass is not fixed but A is, this means that one is in effect considering situations 
in which there are exchanges of energy between the matter surrounding the black hole 
and infinity while leaving the black hole itself unchanged. On the other hand if one fixes 
Madm while letting A vary, one is comparing situations in which the black hole and the 
surrounding matter exchange energy, but no energy is exchanged with infinity. 

Clearly the process of black hole evaporation belongs to the second situation. There- 
fore in the next section we shall consider the action 

S = Jdt {n ab g ab +pq-NH- WH^ + 6A/8tt (22) 

This action gives the time evolution of matter as a function of O for different values of 
A with Madm fixed. 



4 Gravitational Back Reaction to Hawking Radia- 
tion 

In this section we calculate the new expressions of the different wave functions introduced 
in section ^| by replacing the matter action in the given Schwarzschild geometry by the 
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gravitational action eq. (p2|). We then to use these wave functions to compute the 
gravitational corrections to black hole radiation. 

We start with the description of black hole radiation based on the readings of a static 
detector. It is through the change of the deformation of the geometry induced by the 
change of the quantum state of the detector that the change in area will appear in the 
transition rates. 

We first compute the new time dependence of the wave functions associated with the 
two states of the detector. Since the detector is at r = const, the geometry is static and 
the pq and ft ab g a b terms in the action vanish. Moreover, on-shell, the constraints also 
vanish. Hence the only term contributing to the total action (matter + gravity) is the 
surface term at the horizon. This term is equal to QA(g)/8n or QA(e)/8ir where A(g) 
(A(e)) is the horizon area when the detector is in the ground (excited) state. Thus the 
time dependence of the wave functions are 

^(9) = e ieA{9)/8lT 

^ e (6) = e ieA{e)/8n (23) 

It remains to determine the new expression for the outgoing modes which replaces 
eq. (||). As in the absence of back-reaction, the structure of the modes must be such that 
they determine the vacuum near the horizon. Since has been defined at the horizon, 
its relationship to the inertial light like coordinates Uh or , Vhor also defined at the horizon 
is of the form Uh or = — e _e . This encodes the exponential Doppler shift which is the 
hall mark of horizons. Therefore, the new expression is 

0.(6, r) = De lC ^~ B (24) 

in the place of e lCwe * see eqs. (||, §). 

As in Section 2, the transition amplitudes are given by the "time" integral of the 
product of the three waves. Up to the same overall constant, they are given by 



A g ^ e+U) = J dG v s (e)Ve(©)*^(©)* 

A e ^ g+U> = I dG V ff (©)U(©)^(©)* (25) 

Note that the area of the horizon decreases if the detector is in its excited state. There- 
fore, in the limit of small AE, i.e. in the test particle limit, the factor i/j g (Q)*ijj e (Q) = 
e ie(A(e)-A(g))/8n tends ^ Q fa e background field expression e~ ltAE with the correct sign of 
the phase. 

Using the new expressions eq. by replacing by G + in in the upper amplitude, 
one obtains 



\A ' 2 



g^M_ = e -{A(e)-A(g))/A ^6) 



It is thus {A(e) — A(g))/A, the difference of the horizon areas if the detector is in its 
excited or ground state which governs the equilibrium distribution of the detector's 
states. This describes a micro-canonical distribution since we are considering exchanges 
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of energy between the black hole and the detector with no mass exchange at infinity. 
This distribution replaces the canonical expression of eq. (|TlD, governed by the energy 
change E e — E g and Hawking's temperature n/2n. 

The new ratio of the transition rates is a function of the change in area because the 
integrand of the transition amplitude contained the product ifj*ip e which describes the 
replacement of one classical solution by the other one. It is indeed through this 
product that the change in area entered into the expression since its phase is given by 
the difference of the total matter+gravity actions. As noted above this difference tends 
to —AEt when AE is small. This recovery of the background field wave functions in the 



limit of small differences of matter energy is a generic feature |22j of taking into account a 
neglected degree of freedom (here gravity). The same mechanism arises upon taking into 
account recoil effects in scattering amplitudes, see 0. It also explains why one recovers 
the conventional expressions of transition amplitudes starting from the solutions of the 
Wheeler-DeWitt equation in quantum cosmology 0. 

In this calculation the change in area is due to the change of the quantum state of 
the detector. However the existence of a detector is not intrinsic to black hole radiation: 
the detector was only used to reveal the existence of the quanta. Therefore we seek for 
an intrinsic derivation of black hole radiance in which the change in area is due to the 
emission process itself. To this end, we must introduce a model for the matter waves 
which takes into account the deformation of the gravitational background. The simplest 
model is that of KKW in which the matter is described as a spherically symmetric (s- 
wave) light-like shell. The main result of || is that outside and inside the shell the metric 
is Schwarzschild, but with the mass parameter Madm and Madm — A respectively where 
A is the energy of the shell. Then the time parameters inside and outside the shell are 
different and related by junction conditions on the shell. In what follows this relationship 
will play no role. The shell follows a light like geodesic in both the inside and outside 
metric. This property allows us to calculate the action of the shell in a straightforward 
manner. 

Inside the shell the metric is 

ds 2 = -(1 - 2M'/r)dt 2 + (1 - 2M'/r)~ 1 dr 2 + dtt 2 (27) 

where M' = Madm — A is the final mass of the black hole. The shell follows a geodesic 
in this metric. 

^ = (l-2M'/r sh ) 

t sh (r,M') = r + 2M'ln(r-2M') = r + — ln(r - 2M') (28) 

2k' 

where k' = 1 /AM'. 

As in the previous treatment, it is particularly convenient to reexpress the evolution 
in terms of the hyperbolic angle 6 using the Jacobian dQ/dt = k': 

e sh (r, M') = k't + ~ ln(r - 2M') (29) 
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Notice that close to the horizon, the description of the shell's trajectory is now expressed 
only in terms of quantities locally defined, i.e. insensitive to the matter distribution at 
larger radii. 

The action of the shell can in principle be calculated by solving the Hamilton- Jacobi 
equations. Since A' is a constant of motion, these can be solved by separation of variables 
to yield the form S(A', 0) = A'Q/8tt + f(r, A'). Furthermore dS/dA' = must give the 
equations of motion which implies that 

df _ Q sh (r,M') 

OA' ~ 8n m 

Inserting eq. (p9|), this can be integrated to yield 

= \ A A ln(r - r hor (A)) + k(A)t) (31) 

J A' 07T I 

where we have introduced rh or {A), the radius of the horizon parameterized by its area. 
Similarly we have introduced k(A). Note however that k is not strictly speaking the 
surface gravity of the black hole, rather it is determined by the form of the metric near 
the horizon eq. (|l9l). 

In order to obtain the action which governs the transition from A to A', we must 
subtract from S(A') the action AQ/8tt which describes the geometry in the absence of 
the shell. This yields 

r A dA 

S shell (A', A) = -G(A - A')/8n + / — ln(r - r hor (A)) + ... (32) 

J A 1 107T 

where we have dropped the non logarithmic term which is unimportant near the horizon. 
(Note that A\A = A — A' is always positive). Upon postulating that the wave functions 
are WKB, the Schwarzschild mode describing the shell is <^a,+ — e lS ^ A when the 
gravitational interaction is taken into account 

Our task is to now calculate the Bogoljubov transformation between Unruh modes 
and Schwarzschild modes. To impose the regularity of the Unruh modes at the horizon 
we use the Eddington-Finkelstein coordinates inside the shell v = t+(j-? ln(r — r hor ) +r) 
and r, and then impose analyticity in the upper complex r plane. In these coordinates 



eq. fl32|) becomes 

S shM = —A\Ak'v/87c + ^ ln(r - 2M') + t ^- ln(r - r hor (A)) (33) 

l07T J A' 1D7T 

up to non logarithmic terms. Analyticity in the lower half plane imposes that the action 
on the right differs from the action on the left by taking r — 2M' — > (r — 2M')e l7r , hence 
ln(r — 2M') ln(r — 2M') + in . Since both the second and the third term equally 
contribute, we obtain 

S left = S right + i{A x A/8) (34) 
The Bogoljubov transformation between Unruh and Schwarzschild modes is therefore 

A = iV(A)(^, + + e- A ^/VA,-) (35) 
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corresponding to the Bogoljubov coefficients 

I/?a| 2 



-A A A/4 



\oc\ 



2 



(36) 



and one recovers eq. (0]). 

Note that the two derivations lead to the same characterization of the vacuum at the 
horizon. There is indeed a linear combination of positive frequency Unruh modes <p\ eq. 
(|35|) which at fixed r has the dependence of the form De lCuje 



5 Discussion 

We have shown that upon taking into account the gravitational back reaction, the prob- 
ability for a black hole to emit a quantum is given by the exponential of the change in 
area of the black hole. The appearance of this factor AA has the same origin as the 
appearance of the term dA in the first law of black hole mechanics, namely the surface 
term QA/8tt at the horizon. This derivation is not only valid for Schwarzschild black 
holes. Indeed both derivations are also valid for charged black holes and should be easily 
generalizable to rotating black holes. Furthermore the first derivation which does not 
make appeal to any symmetry also applies to the transitions of uniformly accelerated 
detectors in Minkowski space (the Unruh effect). 

In obtaining this result we have made an assumption about the physics at Planckian 
scales: we have postulated that the specification of the vacuum state can be implemented 
by the usual condition of regularity at the horizon when expressed in terms of the local 
coordinates 0, Uhor and r. In other words, we have assumed that the back-reaction does 
not destroy the usual analytical characterization of the vacuum. One would hope to have 
a justification for this assumption at a more fundamental level. This probably requires 
a microscopic description of the physics at the horizon at the Planckian scale, either 
string theory or some other quantum theory of gravity. 
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